Measurement of thin films by polarized light

ABSTRACT

In order to measure the thickness and refractive index of a thin film on a substrate, such as a film of silicon dioxide on a substrate of silicon, a beam of substantially monochromatic polarized light is directed on the film. The reflected light is transmitted through an optical compensator and an optical analyzer both of which are rotating at different angular speeds, ω A  and ω C , respectively; and the transmitted optical intensity is measured as a function of time. A Fourier analysis, for example, of the profile of this optical intensity vs. time can then be used for determining the Stokes parameters of the light reflected by the thin film and thereby also the thickness and refractive index of the film.

BACKGROUND OF THE INVENTION

1. Field of the Invention

This invention relates to the field of the measurement of the characteristics of thin films and, more particularly, to the measurement of the thickness and refractive index of thin films by means of optical techniques.

2. Description of the Prior Art

In order to measure the thickness and refractive index of a thin film on a substrate, various optical techniques have been employed in the prior art. Ordinarily these techniques involve measurements of the intensity and polarization of the transmitted or reflected light from the thin film, as a function of angular orientation ("ellipsometry"). However, the use of these techniques generally tends to be unduly time-consuming. While the more rapid techniques using rotating analyzer ellipsometry have been used in the prior art, such techniques suffer from ambiguities and do not yield correct results in the presence of diffuse reflectance.

In U.S. Pat. No. 3,873,209 issued to Schinke et al. on Mar. 25, 1975, an optical thin film waveguide technique was described for measuring the thickness and refractive index of a thin film. In that technique, an evanescent optical wave is coupled into the film; and the resulting optical intensity minima, characterizing the light which is coupled out of the film as a function of angle, yield data from which the desired thin film thickness and refractive index can be determined. It would be desirable, however, to have a relatively fast thickness measurement technique which does not require the operator to locate such minima.

SUMMARY OF THE INVENTION

In order to measure the thickness and refractive index of a thin film on a substrate, a substantially monochromatic beam of linearly polarized light is incident on the film. The beam undergoes multiple reflections both at the interface of the film with the ambient atmosphere and at the interface of the film with the substrate. The invention is applicable to the cases where either, both, or none of these reflections has a diffuse (nonspecular) component. The resulting reflected beam from the thin film-substrate system is then allowed to pass through a rotating optical compensator and an advantageously faster rotating optical analyzer. The optical beam intensity transmitted by the rotating analyzer is measured by an optical detector as a function of time. An analysis of this function of time yields the Stokes parameters of the light reflected by the thin film-substrate system, from which both the thickness and refractive index of the thin film can be determined.

In a specific embodiment of the invention, a beam of substantially linearly polarized monochromatic light is directed onto the surface of a thin film on a substrate. The film need not be completely transparent and there may be multiple diffuse reflections occurring in the film. The reflected optical beam from the thin film-substrate system is passed through, in spatial sequence, a rotating optical compensator and a rotating optical analyzer, the analyzer rotating at substantially three times the angular speed of the compensator. The optical beam intensity transmitted by the rotating analyzer is measured as a function of time. This function is then Fourier analyzed, and the Fourier coefficients are then used to calculate the desired refractive index and thickness of the thin film.

BRIEF DESCRIPTION OF THE DRAWING

This invention together with its features, objects, and advantages may be better understood from the following detailed description when read in conjunction with the drawing in which FIG. 1 illustrates a thin film being tested in accordance with a specific embodiment of the invention; and

FIG. 2 illustrates a thin film being tested in accordance with another specific embodiment of the invention.

DETAILED DESCRIPTION

As shown in FIG. 1, a source 10 of substantially monochromatic optical radiation directs its output radiation through an optical (nonrotating) polarizer 11, to produce a linearly polarized monochromatic optical beam 12. This beam 12 is incident at an angle φ to the normal of an exposed major surface of a thin film-substrate system 20 comprising a thin film 21 on a substrate 22. For example, the thin film is essentially silicon dioxide and the substrate is essentially monocrystalline semiconductive silicon. Upon multiple internal reflections in the film 21, a reflected beam 13 emerges and passes through an optical phase compensator 14. This compensator is designed so as to introduce a relative phase delay δ, typically of about π/2 radians, between a pair of mutually perpendicularly polarized optical beam components. While the beam 13 passes through, the compensator 14 is rotated at an angular velocity ω_(C) about an axis substantially parallel to the direction of propagation of the beam. A beam 15 is thus transmitted by, and emerges from, the compensator 14 and is then incident on an optical analyzer (polarizer) 16 rotating at an angular velocity ω_(A) about the same axis as that of ω_(C). A beam 17 is thus transmitted by the rotating analyzer 16 and is incident upon utilization means 18 for detecting the optical intensity substantially continuously as a function of time. The utilization means thereupon also performs a Fourier analysis of this optical intensity vs. time, and processes the analysis to yield the refractive index and thickness of the thin film 21 as discussed below in greater detail.

As an example, for the purpose of illustration only, the optical source 10 produces substantially monochromatic visible light, for example, a filtered xenon arc yielding typically a 6000 Angstrom beam with a wavelength spread of advantageously less than 10 Angstroms, perferably less than 5 Angstroms. The azimuth angle P of polarizer 11 is oriented so that the plane of the electric vector associated with the linearly polarized beam 12 is at an angle of about 30° with respect to the plane of incidence determined by the propagation direction of the beam 12 and the normal N to the exposed surface of the thin film 21. This thin film is essentially silicon dioxide which has been grown on a silicon substrate so that the beam undergoes multiple reflections at the ambient-oxide interface and at the substrate-oxide interface and emerges as the reflected beam 13. The compensator 14 is a standard quarter-wave plate rotating at angular speed ω_(C) which is one-third the angular speed ω_(A) of the analyzer 16. The analyzer 16 and the compensator are advantageously rotated around a common axis by a pair of drive chains driven by a single motor, the drive chains having relative tooth or gear ratios of 3:1 for this purpose of obtaining ω_(A) = 3ω_(C). Typically, ω_(C) itself is about 150 radians per second, but other speeds are obviously useful in this invention.

In order to understand in general the selection of optimum parameters, or approximations thereto, which are useful in the practice of the invention, the subscript σ will be used to describe those parameters which are associated with the component of the various optical beams which has its electric vector perpendicular (German "senkrecht") to the plane of incidence, whereas the subscript π will be used for those parameters associated with electric vector parallel to the plane of incidence. The optical amplitude reflectances (which may be complex) for the σ and π components of the beam 12 are denoted by r.sub.Σ and r.sub.π, respectively, with respect to the three-phase system of ambient-thin film-substrate. It is convenient to define a polarization reflectance ratio ρ as the ratio of r.sub.π to r.sub.Σ, which can be a complex number:

    ρ = r.sub.π /r.sub.Σ  = tanψe.sup.i       (1)

This reflectance ratio ρ is of much importance in the determination of the refractive index and thickness of the thin film 21, as discussed below.

In the case of a simple two-phase system, with both phases characterized by real number scalar dielectric constants, then the reflectance ratio ρ is a real number, so that the argument Δ of the exponential of ρ (which should be conveniently limited to range 0 to -π) is either 0 or -π radians, changing suddenly from 0 to -π as a function of angle of incidence φ when φ goes through Brewster's angle. More generally (including cases for complex ρ), the principal angle of incidence φ_(P) is defined as that angle at which Δ = -π/2 (or jumps suddenly from 0 to -π), and the pseudo-Brewster angle φ_(B) is defined as that angle of incidence at which tan ψ (the real part of ρ) is a minimum. The difference between φ_(B) and φ_(P) is ordinarily quite small, and these angles are of interest in this invention in that the principal angle φ_(P) corresponds to that value of the angle of incidence of the two-phase substrate-ambient system for which the optimum accuracy can be obtained in the measurement of the refractive index of the thin film 21 on the substrate 22 in the three-phase system of ambient-oxide-substrate. However, this accuracy is not very sensitive to the angle of incidence, φ, so that useful results can be obtained even if φ deviates from the optimum value φ_(P) by as much as 10° or 20° or more.

While the invention can be practiced with a variety of relative values of angular speeds ω_(A) and ω_(C), it is convenient to use ω_(A) = 3ω_(C), that is:

    A = ω.sub.A t = 3ω.sub.C t                     (2)

    C = A/3 + C.sub.O                                          (3)

where C is the azimuth angle of the major (fast polarization) axis of the compensator 14 relative to the plane of incidence, and A is the instantaneous azimuth angle of the polarization component which is transmitted by the analyzer 16 relative to the plane of incidence. Thus, at t=0, the analyzer transmits the polarization component whose electric vector is in the plane of incidence (plane of the Figure itself); whereas the compensator at t=0 is oriented such that the fastest polarization of electric vector is at an angle of C_(O) with respect to the plane of incidence. For ease of operation and calculation, it is convenient but not essential to set C_(O) = 0 or at least not to exceed 30 degrees. With this setting of C_(O), the intensity I (normalized to an arbitrary value I_(O)) of the transmitted beam 17 detected by means 18 will be a functon of time. It is this normalized intensity (I/I_(O)) which, as a function of time contains the desired information to be processed, by Fourier analysis as described below for example, for yielding the desired value of reflectance ratio ρ and thence the refractive index and thickness of the thin film 21.

Theory

The instantaneous amplitude of the optical electric vector E_(t) of the transmitted optical beam 17 is given by known principles: ##EQU1## where E.sub.π and E.sub.Σ are the instantaneous optical electric vectors respectively parallel to and perpendicular to the plane of incidence, and where δ is the known (measured) phase delay associated with the optical electric vector along the minor (slow) axis of the compensator 14 relative to the optical electric vector along the major (fast) axis of the compensator. The normalized intensity (I/I_(O)), which is equal to 2| E_(t) |² /(|E.sub.π |² + | E.sub.Σ |²), is then given by: ##EQU2## where s_(O), s₁, s₂ and s₃ are the so-called Stokes parameters, defined by:

    s.sub.O = E.sub.π *E.sub.π   + E.sub.Σ *E.sub.Σ (6a)

    s.sub.1 = E.sub.π *E.sub.π  - E.sub.Σ *E.sub.Σ (6b)

    s.sub.2 = 2Re E.sub.π   E.sub.Σ *                 (6c)

    s.sub.3 = -2Im E.sub.π   E.sub.Σ *.               (6d)

Substituting in Equation (4) the expression for C given by Equation (3), and defining a convenient quantity θ as:

    θ = 2A/3 = 2ω.sub.C t - C.sub.O ,              (7)

it is found from Equation (5) that: ##EQU3## It is convenient to introduce Fourier coefficients α₁, α₂, α₃, β₁, β₂, β₃, defined by: ##EQU4## The values of all these Fourier coefficients can be measured by the utilization means 18, keeping in mind the definition of θ as a function of time given by equation (7). Comparing equations (8) and (9), it is found that the relative Stokes parameters are given by:

    s.sub.1 /s.sub.O = α.sub.3 [ 1+(α.sub.1.sup.2 +β.sub.1.sup.2).sup.1/2 /(α.sub.3.sup. 2+ β.sub.3.sup. 2).sup.1/2 ]                                              (10a)

    s.sub.2 /s.sub.O = β.sub.3 [1+(α.sub.1.sup.2 +β.sub.1.sup.2).sup.1/2 /(α.sub.3.sup.2 +β.sub.3.sup.2).sup.1/2 ]                            (10b)

    s.sub.3 /s.sub.O = -(α.sub.2.sup.2 +β.sub.2.sup.2).sup.1/2 /sinδ sgn(β.sub.2)                             (10c)

where only positive square roots (1/2 powers) are taken. Thus, the Fourier analysis of I/I_(O) yields the values of the relative Stokes parameters. From these the desired value of reflectance ratio ρ can be calculated from: ##EQU5## where P is the actual (measured) value of the azimuth angle of the polarized 11, and i is the imaginary root. It should be noted at this point that the prior art ellipsometry techniques involving rotating only the analyzer assume that (s₁ /s_(O))² +(s₂ /s_(O))² +(s₃ /s_(O))² is equal to unity, which is not true in the presence of diffuse reflection. In addition, these prior art techniques determine only the magnitude of s₃ /s_(O) but not its algebraic sign (s_(O) being always positive: Eq. (6a)), thereby producing the ambiguity of two possible values of the polarization reflectance ratio ρ instead of solely the correct one. Thus, this invention yields unambiguously the single value of reflectance ratio ρ. From this value of ρ, the refractive index and thickness of the thin film 21 can be calculated, as set forth in a paper by A. R. Reinberg entitled "Ellipsometer Data Analysis with a Small Programmable Desk Calculator" in Applied Optics, Vol. 11, No. 5, pp. 1273-1274 (May 1972). Briefly, the calculation involves a method of successive approximation (Newton's method), starting with the known relationship of the reflectance ratio ρ to the reflection parameters of the thin film-substrate system: ##EQU6## where, for example, r.sub.π ,_(fs) is the reflection coefficient of the π polarization of the two-phase film-substrate interface, the subscript pair af refers to the two-phase ambient-film system, and where: ##EQU7##

    Z = exp[4πi(ε.sub.f -ε.sub.a sin.sup.2 φ).sup.1/2 d/λ]                                               (17)

where n_(a) is the known refractive index of the ambient, φ is the actual angle of incidence, and where ε_(f) is the unknown thin film dielectric constant relative to vacuum, d is the unknown thickness of the film, ε_(s) is the known substrate dielectric constant relative to vacuum, ε_(a) =n_(a) ² is the known ambient dielectric constant, and λ is the known vacuum wavelength of the substantially monochromatic light. Accordingly, there are two ultimate unknowns, ε_(f) and d, and the intermediate unknown Z. From equation (17), however, it is seen that Z must have an absolute magnitude of unity (|Z|=1) (so long as ε_(f) =n_(f) ² is greater than ε_(a) sin² φ = n_(a) ² sin² φ, which will always be true for usual cases of interest involving a thin film whose refractive index is greater than that of the ambient). Moreover, putting into equation (12) the measured value of reflectance ratio ρ, and substituting into equation (12) the various expressions for the reflection coefficients of equations (13-16), then a quadratic equation in Z results containing only ε_(f) as the only other unknown. In order to solve this equation, a first trial value of ε_(f) is assumed and a consequent first approximation to both Z and d|Z| /dε_(f) is calculated for this trial value of ε_(f), the first approximation to Z being selected as the root of the quadratic equation in Z which is closer to |Z|=1, and d|Z|/dε_(f) being determined therefrom by known differential calculus methods. Then a second trial value of ε_(f) is adopted using the straight line slope approximation of d|Z| /dε_(f) and the above first approximation of Z (Newton's method) to try to make | Z|=1 with this second trial value. A third trial value of ε_(f) is then adopted if the magnitude of the second approximation of Z (root closer to | Z| =1) on the second trial is not sufficiently close to unity within the desired accuracy, then using the second value of ε_(f) for the d|Z|/dε_(f) straight line approximation. For each (complex) value of Z thus found, using the corresponding trial value of ε_(f), then d can be found from equation (17). When the trial values ε_(f) and the corresponding values of d converge sufficiently to a limit, in accordance with desired accuracy, then the process of successive approximation is terminated, and the last trial value of ε_(f) = n_(f) ² provides the desired refractive index, n_(f), of the thin film, whereas the value of d corresponding to this last trial value of ε_(f) is the desired thickness of the thin film.

While the invention has been described in detail in terms of a specific embodiment, various modifications can be made without departing from its scope. The substrate can be any material which reflects light from the source beam 12, whether the reflection is specular, diffuse, polarizing, or depolarizing, such as glass, metal, plastics or liquids. The thin film can be any material which forms a transparent or semitransparent layer on the substrte. Determination of the thickness and dielectric constant of films which are not transparent is also possible but requires the use of measurements, at more than a single wavelength or angle of incidence. The invention is not restricted to angular speeds ω_(A) and ω_(C) with ω_(A) /ω_(C) =3, but other relative speeds can be used such as ω_(A) /ω_(C) =4, 5, 6, . . . , or 7/2, 9/2, 11/2, . . . , or other rational number; but it should be noted that certain ratios or ω_(A) /ω_(C), such as ω_(A) /ω_(C) = 2, should be avoided, in order to reduce undesired ambiguities.

In the case where there is no diffuse reflectance, the compensator can be placed before (instead of after) the reflection, along the propagation direction, in the practice of the invention (as indicated in FIG. 2). Moreover, instead of linear polarization, the polarizer 11 can produce other known types of polarization, such as circular or elliptical.

It should be understood that there is a single remaining ambiguity in the calculation of the thickness d from Z in equation (17) due to the periodicity of the imaginary exponential, thereby yielding a spectrum of discrete possible values of d. Such an ambiguity can be removed, for example, by using a second different wavelength and obtaining a second discrete spectrum of values of d and selecting the value of d common to the two spectra. 

What is claimed is:
 1. A method of testing a thin film-substrate system which comprises the steps ofa. directing an incident beam of substantially monochromatic polarized optical radiation onto an exposed surface of the thin film, and b. detecting the resulting reflected beam of radiation reflected from the surface and transmitted through a rotating optical compensator and a rotating optical analyzer, the compensator and analyzer rotating about an axis along the propagation direction of the reflected beam at different first and second different angular speeds, respectively.
 2. The method of claim 1 in which the step of detecting includes measuring the intensity of resulting beam as a function of time and in which the incident beam is substantially linearly polarized.
 3. The method of claim 2 in which the reflected beam passes through the rotating compensator prior to passing through the rotating analyzer.
 4. The method of claim 3 wherein said compensator is rotating at a lower angular speed than that of the analyzer.
 5. The method of claim 4 in which the compensator is rotating at substantially one-third the angular speed of the analyzer.
 6. A method of testing a body including a thin film which comprises the steps ofa. directing a substantially linearly polarized monochromatic beam of optical radiation onto an exposed surface of the thin film whereby a reflected beam is produced by reflection at said surface; b. rotating an optical compensator in the path of the reflected beam at a predetermined first angular speed about an axis parallel to the propagation direction of the reflected beam, while said reflected beam traverses through said compensator; c. rotating an optical analyzer in the path of the beam coming from the rotating compensator at a predetermined second angular speed different from the first angular speed about the said axis; and d. detecting the intensity of the beam transmitted by the rotating analyzer as a function of time.
 7. The method of claim 6 in which said detecting is performed over a time interval including at least one complete rotation of the compensator.
 8. The method of claim 6 in which the angular speed of rotating the analyzer is greater than that of the compensator.
 9. The method of claim 8 in which the angular speed of rotating the analyzer is substantially three times that of the compensator.
 10. A method of testing a thin film on a substrate body which comprises the steps of:a. directing a beam of polarized substantially monochromatic optical radiation through a rotating optical compensator onto a surface of the thin film; and b. detecting the optical radiation reflected by the body after the reflected radiation passes through an optical analyzer rotating at an angular speed which is substantially different from the angular speed of the rotating compensator.
 11. The method of claim 10 in which the reflected radiation is detected substantially continuously for a period of at least a full rotation of the analyzer.
 12. The method of claim 11 in which the beam incident on the rotating compensator is substantially linearly polarized.
 13. The method of claim 12 in which the analyzer is rotating at substantially three times the angular speed of the compensator. 